Is Earthquake Triggering Driven by Small Earthquakes? 
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Using a catalog of seismicity for Southern California, we measure how the number of triggered 
earthquakes increases with the earthquake magnitude. The trade-off between this relation and 
the distribution of earthquake magnitudes controls the relative role of small compared to large 
earthquakes. We show that seismicity triggering is driven by the smallest earthquakes, which trigger 
fewer events than larger earthquakes, but which are much more numerous. We propose that the non- 
trivial scaling of the number of triggered earthquakes emerges from the fractal spatial distribution 
of seismicity. 
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Large shallow earthquakes are always followed by af- 
tershocks, that are due the stress change of the main- 
shock. The number n{M) of aftershocks of a mainshock 
of magnitude M has been proposed to scale with AI as 
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This relation accounts for the fact that large earthquakes 
trigger many more aftershocks than small earthquakes. A 
similar relation holds for the distribution of earthquake 
magnitudes P{M) ^3 given by 
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with b typically close to 1, which implies that small earth- 
quakes are much more frequent than large earthquakes. 

Because large earthquakes release more energy and 
trigger more aftershocks than smaller earthquakes, it is 
usually accepted that interactions between earthquakes 
and earthquake triggering are dominated by the largest 
earthquakes. However, because they are much more fre- 
quent that larger earthquakes, small earthquakes are also 
just as important as large earthquakes in redistributing 
the tectonic forces if 6 = 1 0|. Other quantities, such as 
the Benioff strain e lO"'^^*^, are dominated by small 
earthquakes if 6 > 0.75. 

The a-exponent is an important parameter of earth- 
quake interaction that is used in many stochastic models 
of seismicity or prediction algorithms 0, IE S S S Qi 
Il7j . This parameter controls the relative role of small 
compared to large earthquakes. While there is a signifi- 
cant amount of literature on the 6- value, very few studies 
have measured accurately the a exponent in real seismic- 
ity data. M any studies use a — h without justification 
0, la llJ; ll3 • In this case, small earthquakes are just 
as important as larger ones for the triggering process. Us- 
ing ^ and the global number N{M) of aftershocks 
triggered by all earthquakes of magnitude M scales as 



N{M) = P(M) n(M) - 10^"-^) 
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and is indeed independent of M in the case a = b. In 
the case a < b, aftershock triggering is controlled by 



the smallest earthquakes, while the largest earthquakes 
dominate if a > 6. 

A few studies measured directly a from aftershocks se- 
quences, using a fit of the total number of aftershocks 
as a function of the mainshock magnitude 0, 0, UlL Il3j | . 
These studies yield a- value close to 1, but the limited 
range of the mainshock magnitude considered and the 
large scatter of the number of aftershocks per mainshock 
do not allow an accurate estimation of a. The case 
a — b also explains another well documented property 
of aftershocks, known as Bath's Law 0,0,^8^, which 
states that the difference between the mainshock mag- 
nitude and its largest aftershock is on average close to 
1.2, independently of the mainshock magnitude. Again, 
the limited range of mainshock magnitudes used in these 
studies and possible biases of data selection [l^ does not 
allow to test the dependence of the magnitude difference 
as a function of the mainshock magnitude. 

Other studies measured a indirectly using a stochas- 
tic triggering model called "Epidemic Type Aftershock 
Sequence" model (ETAS) H 0, 113 based only on the 
Gutenberg-Richter and Omori laws 0, 0|- This model 
assumes that each earthquake above a magnitude thresh- 
old mo can trigger direct aftershocks, with a rate that in- 
creases as ~ 10"*^ with its l uag nitude, and decays with 
time according to Omori law j21| . The average total num- 
ber of aftershocks n(M), including the cascades of indi- 
rect aftershocks, has the same dependence ~ iQ"*^ with 
the mainshock magnitude M as the number of direct af- 
tershocks. Using this model, a can be measured using a 
maximum likelihood method 0,0 ESl- For instance, [iOl 
analyzed 34 aftershock sequences in Japan and measured 
a in the range [0.2 — 1.9] with a mean value of 0.86. The 
a-values obtained from the inversion of this model are 
not well constrained due to the small number of events 
available and to possible biases of the inversion method. 
Indeed, these studies do not take into account the influ- 
ence of earthquakes below the detection threshold, which 
may significantly bias the estimation of a. This method 
may also be biased by the incompleteness of the catalog 
just after the mainshock, and by possible trade-offs be- 
tween the ETAS parameters. The regime a > 6 of the 
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ETAS model is probably not relevant for real seismicity. 
If we do not assume a roll-off of the magnitude distribu- 
tion P{M) for large M , this regime gives a finite time 
singularity of the seismicity rate which goes to infinity 
in finite time as l/(tc — ^)'" H^- Such a power-law 
increase of seismic activity can describe the acceleration 
of the deformation preceding material failure as well as a 
starquake sequence 23] , but cannot describe a stationary 
seismic activity. 

In this study, we use a stacking method to estimate 
the average rate of earthquakes triggered (directly or in- 
directly) by a previous earthquake as a function of the 
magnitude of the triggering earthquake. We use the seis- 
micity catalog of Southern California provided by the 
Southern California Data Center [2^, which covers the 
time period 1975-2003, and which is complete above 
M = 3 for this time period. The magnitude distribu- 
tion shown in Fig. lb follows the Gutenberg-Richter 
law for M > 3 and attests for the completeness of the 
catalog for this time period above magnitude 3. We 
do not use the usual distinction between "foreshocks" , 
"mainshocks" and "aftershocks" , and the constraint that 
the aftershocks must be smaller than the mainshock be- 
cause this classification is not based on physical differ- 
ences. Indeed, recent studies have shown that a sim- 
ple model that assumes that each earthquake can trigger 
earthquakes of any magnitude, without any distinction 
between "foreshocks", "mainshocks" and "aftershocks" 
can reproduce many properties of real seismicity includ- 
ing realistic foreshock sequences @,H0,0,|23,|23- Con- 
straining "triggered earthquakes" to be smaller than the 
mainshock would obviously underestimate the number 
of earthquakes triggered by small mainshocks and thus 
overestimate a. 

We define a "triggered earthquake" as any event oc- 
curring in a space-time window R x T after a preceding 
"mainshock" above the background level, whatever the 
relative magnitude of the triggered and triggering earth- 
quakes. We consider as a potential "mainshock" each 
earthquake that has not been preceded by a previous 
larger earthquake in a space-time window d x T in order 
to estimate the rate of seismicity triggered by this "main- 
shock" removing the influence of previous earthquakes. 

This definition of "triggered earthquakes" and "main- 
shocks" contains unavoidably a degree of arbitrariness in 
the choice of the space-time windows but the estimation 
of a is found to be robust when changing T, R and d. 
We have tested different methods for the choice of R, ei- 
ther fixed or increasing with the mainshock magnitude. 
We use a distance R increasing with the mainshock mag- 
nitude because the aftershock zone is usually found to 
scale with the rupture length '2^ . We use R equal to 
1 rupture length of the mainshock. For small mainshock 
magnitudes, this choice would lead to unacceptable val- 
ues of R smaller than the location error, and thus to un- 
derestimate the number of triggered earthquakes of small 



mainshocks. Therefore, we impose R > 5 km, larger than 
the location error. Taking R fixed has the advantage of 
not introducing by hand any scaling between the after- 
shock zone and the mainshock magnitude. However, it 
may overestimate the number of earthquakes triggered 
by the smallest mainshocks if R is too large, or under- 
estimate the number of triggered events of the largest 
mainshock if R is too small. 

The results obtained for T = 1 year, R = 0.01 x IQ^-^^i 
km and c? = 50 km are presented in Fig. 1. The rate 
of triggered earthquakes is found to decay according to 
Omori's law K(M)/tP, with the same exponent p w 0.9 
for all mainshock magnitudes M (Fig. la). The ampli- 
tude K{M) increases exponentially ~ IQ"*'^ as a function 
of M with a = 0.81 (Fig. lb). This confirms that the 
scaling of the rate of triggered earthquakes with M fol- 
lows (^. Our method is more accurate than previous 
studies nil HIES- Indeed, these studies 011 El El 
determine the scaling of n(M) with M using the total 
number of aftershocks 0, H ITlL 111 in a time window 
[0 — T] after a mainshock, and can thus be biased by the 
incompleteness of the catalogs just after a large main- 
shock or by the background seismicity at large times af- 
ter a mainshock. In contrast, in order to deal with these 
problems, we use the seismicity rate in the time window 
where we observe the Omori law decay characteristic of 
triggered seismicity. 

The value of a is robust when increasing or decreasing 
the distance R used for the selection of triggered earth- 
quakes between 1 to 5 rupture lengths, or when increas- 
ing the minimum value of R from 2 to 10 km. Selecting 
earthquakes within a disk of fixed radius i? = 50 km 
for all mainshock magnitudes yields a slightly smaller 
value a = 0.72. Decreasing R leads to a smaller value of 
a because it underestimates the number of events trig- 
gered by the largest mainshocks, which have a rupture 
size larger than R. When increasing R from 10 to 100 
km, the value of a first increases with R and then sat- 
urates around a = 0.72 for i? > 50 km. We have also 
checked that a is not sensitive to the parameter d used 
for the selection of mainshocks if d > 50 km. All values 
of a, for reasonable values of the parameters R in the 
range 30 — 100 km, T between 0.1 and 2 yrs and d > 50 
km, and for different time periods of the catalog, are in 
the range [0.7 — 0.9]. We have also tested the method on 
synthetic catalogs generated with the ETAS model. We 
recover the a parameter with an error smaller than 0.05. 

For the same catalog of seismicity, we measure us- 
ing a maximum likelihood method the 6-value of the 
Gutenberg-Richter law © equal to & = 1.08 ± 0.10. We 
have also tested that the magnitude distribution P{M) 
of triggered events is independent of the mainshock mag- 
nitude. Fig. 2 shows the magnitude distribution of trig- 
gered events for different ranges of the mainshock magni- 
tude, using the same data as in Fig. 1. This figure shows 
that a large earthquake can be triggered by a smaller 
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earthquake. Our results suggest that a is significantly 
smaller than the h exponent of the magnitude distribu- 
tion. Whether the exponent a varies with region and 
maybe even with time is an interesting question that is 
outside the scope of this Letter but we urge further stud- 
ies in that direction. 

We now propose a simple explanation for this non- 
trivial scaling of the number of triggered earthquakes 
with the mainshock magnitude, and we suggest that a 
can be related to the fractal structure of the spatial dis- 
tribution of seismicity. It is widely accepted that the 
aftershock zone scales with the rupture length [iil l27j . 
While the area affected by the stress variation induced 
by an earthquake increases with the rupture length, the 
stress drop is independent of the mainshock magnitude 
[islE^ . The stress variation at a distance from the main- 
shock proportional to the fault length L is thus indepen- 
dent of the mainshock magnitude, neglecting the effect 
of the finite width of the crust and the visco-elastic de- 
formation in the lower-crust. Therefore, assuming that 
earthquakes triggered by the stress change induced by 
the mainshock, the density of earthquakes triggered at a 
distance up to i? ~ L from the mainshock is independent 
of the mainshock magnitude. The increase of the number 
of triggered events with the mainshock magnitude results 
only from the increase in the aftershock zone size with 
the rupture length. 

The rupture length is usually related to the magnitude 
by Hi 
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The same relation thus holds between the aftershock zone 
size i? and the mainshock magnitude. 

In order to estimate the scaling of the number of trig- 
gered events with the rupture length, we need to make an 
assumption about the spatial distribution of earthquakes 
around the mainshock. Assuming that triggered earth- 
quakes are uniformly distributed on the fault plane, and 
using 10} , the number of earthquakes triggered by a main- 
shock of magnitude M is given by n{M) ^ 10*^ 9] 
and thus leads to a = 1. The value a — 0.5 obtained 
for a numerical model of seismicity jl2] suggests that in 
this model earthquakes are triggered m ostly on the edge 
of the fracture area of the mainshock Il2|. Our result 
a = 0.8 for Southern California seismicity implies that 
triggered earthquakes are distributed neither uniformly 
on the rupture plane nor on the edge of the rupture, but 
rather on a fractal structure of dimension D < 2. Us- 
ing the definition of the capacity fractal dimension, the 
number of aftershocks is 



n{M) - R 



D 



(5) 



where R is the characteristic length of the aftershock 
zone. Using Q and (0) , we obtain the scaling of the num- 
ber of triggered earthquakes with the mainshock magni- 
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which gives a = 0.5D. Our estimation a = 0.8 for South- 
ern California seismicity thus suggests D — 1.6. This 
value of the fractal dimension of aftershocks hypocenters 
has never been measured for Southern California seismic- 
ity. Our estimate of D is significantly smaller than the 
value measured in the range [2 — 2.8] for aftershock se- 
quences in Japan '20'| . This fractal dimension of the spa- 
tial distribution of triggered earthquakes results in part 
from the fractal structure of the fault system j^O] , but it 
may also reflect the non-uniformity of the distribution of 
the earthquakes on the fault due to the heterogeneity of 
stress or strength on the fault. The fractal dimension of 
the aftershock distribution may thus be smaller than the 
fractal distribution of the fault system. 

While the energy release and the total slip on faults is 
controlled by the largest earthquakes, the suggestion that 
a < b implies that small earthquakes may be more im- 
portant than large earthquakes in triggering earthquakes. 
We have also checked that the magnitude distribution of 
triggered earthquakes is independent of the mainshock 
magnitude (Fig. 2). This implies that earthquake trig- 
gering is driven by the smallest earthquakes at all sc ales, 
even for the largest earthquakes. Other observations |23 
support the conclusion that the same mechanisms can 
explain the triggering of a large earthquake by a smaller 
one and the triggering of a small earthquake by a previ- 
ous larger event. 

Recent studies 0| have proposed that secondary af- 
tershocks dominate an aftershock sequence, so that sub- 
sequent large aftershocks are more likely to be triggered 
indirectly by a previous aftershock of the mainshock. Our 
study further suggests that the smallest earthquakes will 
dominate the triggering of following earthquakes. The 
importance of small earthquakes casts doubts on the rel- 
evance of calculations of direct stress transfer functions 
to predict seismicity |21;], because large earthquakes are 
likely to be triggered by the smallest earthquakes below 
the detection threshold of the seismic network. Small 
earthquakes taken individually have a very low probabil- 
ity of triggering a large earthquake. But because they 
are much more numerous than larger earthquakes, col- 
lectively, they trigger more earthquakes. This result re- 
quires the existence of a small magnitude cut-off mo, 
below which earthquakes may occur but cannot trigger 
earthquakes larger than mg, or a change of the scaling of 
N{M) given by Q for small earthquakes, otherwise the 
seismicity at all scales would be controlled by infinitely 
small earthquakes. 
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FIG. 1: Average rate of triggered earthquakes n{M,t) as a 
function of time t after the triggering earthquake (a) for dif- 
ferent values of the magnitude M of the triggering earthquake 
increasing from 3 to 7 with a step of 0.5 from bottom to top. 
The rate of aftershocks K{M) as a function of M is shown 
in panel (b) (circles) with the cumulative magnitude distribu- 
tion (crosses). K{M) is obtained by fitting each curve n(M, t) 
by K{M)/t°-^ in the range 0.01 < t < 365 days for M < 6.5. 
For large M > 6.5 mainshocks, there is a roll-off of the seis- 
micity rate for small times after the mainshock due to the 
incompleteness of the catalog after large mainshocks, caused 
by the saturation of the seismic network. Therefore we mea- 
sure K{M) in the range t > 0.1 day for M = 6.5 and t > 0.3 
day for M = 7. 
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FIG. 2: Cumulative magnitude distribution of triggered 
earthquakes for different values of the mainshock magnitude 
between 3 (dark line, small circles) and 7 (gray line, large 
symbols) using the same time interval for the selection of af- 
tershocks as in Fig. 



